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We investigate the electronic transport properties of a bilayer graphene flake contacted by two 
monolayer nanoribbons. Such a finite-size bilayer flake can be built by overlapping two semi-infinite 
ribbons or by depositing a monolayer flake onto an infinite nanoribbon. These two structures 
have a complementary behavior, that we study and analyze by means of a tight-binding method 
and a continuum Dirac model. We have found that for certain energy ranges and geometries, 
the conductance of these systems oscillates markedly between zero and the maximum value of 
the conductance, allowing for the design of electromechanical switches. Our understanding of the 
electronic transmission through bilayer flakes may provide a way to measure the interlayer hopping 
in bilayer graphene. 

PACS numbers: 73.22.Pr, 73.23.Ad 



I. INTRODUCTION 

Graphene is a sheet of carbon atoms that order in a 
honeycomb structure, which is composed of two inequiv- 
alent triangular sublattices A and B. Since its exper- 
imental isolation in 2004^ and the subsequent verifica- 
tion of its exotic properties, the interest in this mate- 
rial has boosted. Carriers in monolayer graphene behave 
as two-dimensional (2D) massless Dirac fermionsj 2 ^ with 
a linear dispersion relation e(k) = ±vpk. Phenom ena 
of fundamental nature, such as quantum Hall effeclP^ 
and KleirP^ tunneling have been recently measured in 
graphene based devices. 

Being a material of atomic thickness, graphene is 
regarded as a promising candidate for nanoelectronic 
applications. 2 By patterning graphene, its electronic 
structure can be altered in a dramatic fashion: size quan- 
tization yields ribbons with electronic gaps, essential for 
electronics.^^ By imposing appropriate boundary con- 
ditions, the physics of graphene nanoribbons is well de- 
scribed within a continuum Dirac modelP^iHIU Further- 
more, connections and devices can be designed in a pla- 
nar geometry by cutting graphene layers!^ Another way 
to modify the band structure of graphene is to stack 
two graphene monolayers, 1 and 2, forming a bilayer 
grapheneP^HII] i n bilayer graphene there are four atoms 
per unit cell, with inequivalent sites Al, Bl and A2, B2 
in the first and second graphene layers, respectively. 

Different stacking orders can occur in bilayer graphene. 
Due to its larger stability for bulk graphite, the most 
commonly studied is AB (Bernal) stacking. In the AB 
stacking, the two graphene layers are arranged in such a 
way that the Al sublattice is exactly on top of the sublat- 
tice B2. In the simple hexagonal or A A stacking, both 
sublattices of sheet 1, Al and Bl, are located directly 
on top of the two sublattices A2 and B2 of sheet 2. Al- 
though graphite with direct or AA stacking has not been 
observed in natural graphite, it has been produced by 
folding graphite layers at the edges of a cleaved sample 



with a scanning tunneling microscope tip additionally, 
the growth of AA-stacked graphite on (111) diamond has 
also been reported.^ Furthermore, it has been recently 
found that AA stacking is surprisingly frequent in bi- 
layer graphenej^l so it should be also considered as a 
realistic possibility in few-layer graphene. The interpla- 
nar spacing for the AB stacking has been experimentally 
determined to be cab — 3.35Ap^, whereas for the AA 
stacking seems to be somewhat larger, caa ~ 3.55AP^ 
First-principles calculations agree with these valuespAH^S] 
In any case, the distance between atoms belonging to 
different layers in both stackings is much larger that the 
separation between atoms in the same layer, acc = 1-42 

A. 

Nanostructures based on bilaye r graphene have begun 
to be explored only recently? 16 * 24 ^ ? Bilayer graphene 
nanoribbons might present better signal-to-noise ra- 
tio in transport experiments than monolayer ribbons.^ 
Graphene flakes are quantum-dot-like structures, and be- 
cause of their aspect ratio they are also called nanobars. 
Both, bilayer nanoribbons and bilayer flakes, show in- 
teresting properties with an intriguing dependence on 
stacking. The dependence of the energy gap of bilayer 
graphene flakes on their width and length as well as on 
their atomic termination has been recently reported! 2 ^ 

In this paper we concentrate in the transport proper- 
ties of bilayer armchair graphene flakes with nanoribbon 
contacts. We consider that the most likely way of achiev- 
ing such quasi-zero dimensional structures is either by the 
overlap of two nanoribbons, depicted in the lower part of 
Fig. [T] or the deposition of a finite-size graphene flake 
over a graphene nanoribbon, shown in the upper part of 
Fig. [I] These two configurations correspond to two dif- 
ferent ways of providing monolayer nanoribbon leads for 
the bilayer flake: either the ribbon leads are contacted to 
different layers of the flake, or to the same monolayer. We 
will address these two configurations as bottom-bottom 
(1 — > 1) or bottom-top (1 — > 2), respectively. In both 
geometries the width of the bilayer flake and nanorib- 
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bons is the same, W, and the length of the bilayer region 
is L. In this work we consider narrow armchair metallic 
graphene nanoribbons in the energy range for which only 
one incident electronic channel is active. 



L 




FIG. 1. (color online) Schematic view of two possible geome- 
tries for a bilayer graphene flake contacted by two nanorib- 
bons. Top: A finite-size bilayer graphene flake achieved by 
overlaying a monolayer graphene quantum dot over an infi- 
nite graphene nanoribbon (1 — > 1 configuration). Bottom: 
The bilayer graphene flake is formed by the overlap of two 
semi-infinite nanoribbons (1 — > 2 configuration). In both 
cases the width and length of the bilayer region are L and 
W respectively. 

We calculate the conductance with two different ap- 
proaches: a tight-binding model using the Landauer- 
Biittiker formalism and a mode-matching calculation in 
the continuum Dirac-like Hamiltonian approximation. 
Our main results are the following: 

i) In the AA stacking configuration, the transmission 
through the system shows antiresonances due to the in- 
terference of the two propagating electronic channels in 
the bilayer flake. For a bilayer region of length L we ob- 
tain that the conductance oscillates as function of energy 
with a main period vf^/L. For a fixed incident energy E, 
the conductance as a function of the length L oscillates 
with two main periods: wvf/ji and ttvf/E, being 71 the 
interlayer hopping parameter. The bonding/antibonding 
character of the bilayer bands in the AA stacking makes 
the bottom-top and bottom-bottom conductances to be 
rather complementary: the conductance is zero in the 
bottom-top configuration and it is finite in the bottom- 
bottom arrangement at zero energy, and in general the 
maxima of the bottom-top configuration coincide with 
the minima of the bottom-bottom one and viceversa. 

ii) For the AB stacking, and for energies larger than the 
interplane hopping 71, these devices behave similar to 
those in the AA configuration because there are also two 
propagating electron channels in the bilayer flake at these 
energies. The conductance presents antiresonances with 
periods depending on E, 71, and L. An interesting differ- 
ence is that for a fixed incident energy, the period related 
with interlayer hopping is twice than that found for the 
AA stacking. This reflects that in the AB stacking only 
half of the atoms are connected by interlayer hopping, 



whereas in the A A arrangement all atoms are connected, 
iii) For energies smaller than the interplane hopping, for 
which the AB stacking has only a propagating chan- 
nel, the conductance shows Fabry-Perot-like resonances. 
These are associated with constructive interferences in 
the only available electronic channel. At zero energy the 
conductance of the bottom-bottom configuration is unity, 
whereas in the bottom-top geometry the conductance is 
zero. 

We have analyzed the dependence of the transmission 
with the structural parameters and the interlayer cou- 
pling in bilayer graphene. This study provides a way to 
determine the interlayer hopping by studying the vari- 
ation of the low energy conductance of two overlapping 
nanoribbons with the bilayer flake length; in addition, it 
could clarify the role of stacking in the transport charac- 
teristics of these systems. Our results also indicate that 
the conductance, as function of energy and system size, 
oscillates markedly between zero and a finite value, allow- 
ing for the design of electromechanical switches based on 
overlapping nanoribbons. The introduction of an exter- 
nal gate voltage is of interest for potential applications, 
however, we restrict ourselves to zero gate voltage in or- 
der to obtain analytical expressions in the Dirac model 
and acquire a physical understanding of the transport 
properties of these structures. 

This work is organized as follows. In Section [II] we in- 
troduce the tight-binding and Dirac Hamiltonians we use 
to model the electronic properties of graphene. Section 
|III| is dedicated to describe the conductance calculations, 
both in the tight-binding approximation, for which we 
use Landauer-Buttikcr formalism, and in the continuum 
Dirac-like model, where we use a wavefunction matching 
technique. Section IV is dedicated to present numeri- 
cal results obtained in the tight-binding Hamiltonian and 
compare them with the analytical results obtained in the 
Dirac formalism. Finally, we conclude in Section V sum- 
marizing our main results. 



II. THEORETICAL DESCRIPTION OF THE 
SYSTEM 

The low energy properties in graphene are mainly 
determined by the p z orbitals. Thus, we adopt a tt- 
band tight-binding Hamiltonian with nearest-neighbor 
in-plane interaction given by the hopping parameter 
70 = 2.66 eV. In undoped graphene, the conduction and 
valence bands touch at two inequivalent points of the 
Brillouin zone K and K'. Near these points, the elec- 
tric properties of graphene can be described by a mass- 
less Dirac HamiltoniarP that has a linear dispersion with 

slope VF=^-Jo a o, where do = 2.46 A is the graphene 
in-plane lattice parameter. 

Bilayer graphene consists of two graphene layer cou- 
pled by tunneling. The interlayer coupling is modeled 
with a single hopping 7! connecting atoms directly on top 
of each other, which we take as 71 = O.I70, in agreement 
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with experimental results J2U2QI As discussed in the Intro- 
duction, the interlayer hopping is considerably smaller 
than the intralayer hopping because the nearest-neighbor 
distance between carbon atoms is much smaller than the 
interlayer separation. We do not include other remote 
terms, such as trigonal warping 73, because even though 
it has a similar value to 71 , its effects are more important 
away from the neutrality point, where the Dirac cones are 
distorted and therefore the continuum approximation is 
not so good. 

A. Tight-binding Hamiltonians 

The tight-binding Hamiltonian for the ^IB-stacked bi- 
layer reads 

H AB =- 7o £ (a+^ + h.c) 

<i.j>,rn 

-■yi^2(4,ib2,i + h.c.), (1) 

i 

where a m ,i{bm,i) annihilates an electron on sublattice 
A(B), in plane m = 1,2, at lattice site i. The subscript 
< i,j > represents a pair of in-plane nearest neighbors. 
For the AB stacking we assume that the atoms on the 
A sublattice of the bottom layer (Al) are connected to 
those on the B sublattice of the top layer (B2). The 
second term in Eq. (JlJ represents the hopping between 
these two sets of atoms. 

For the bilayer with AA stacking, all the atoms of layer 
1 are on top of the equivalent atoms of layer 2; thus, the 
Hamiltonian takes the form 

H AA =~lv «A"j ■+ h.c.) 

<i,j>,m 

~7iX!K« a 2.i + 6 M 6 2, l +h.c). (2) 

i 

As we are interested in the transport properties of the 
bilayer flakes, we will concentrate on structures where 
the leads are monolayer armchair graphene nanoribbons 
(aGNR), with widths chosen to have metallic character. 
We denote the ribbon width with an integer N indicat- 
ing the number of carbon dimers along it. With this 
convention, a nanoribbon of width N = 3p + 2, where 
p = 0, 1, 2..., is metallic. In Fig. [2ja) we plot the atomic 
geometry of the monolayer aGNR leads and the corre- 
sponding low energy electronic bands, as obtained from 
the tight-binding Hamiltonian. Note that in aGNR the 
two Dirac points collapse in just onePS Near the Dirac 
point the dispersion is linear, vpk. In the transport cal- 
culations we will only consider incident electrons inside 
this subband, i.e., with energy lower than the second sub- 
band. An aGNR is metallic because of a particular com- 
bination of the wavefunctions coming form the two orig- 
inal Dirac points. This combination is preserved when 
piling up two metallic armchair monolayer ribbons, be- 
ing the corresponding bilayer nanoribbon also metallic. 



The details of the low energy spectrum of bilayer 
nanoribbons depend on the particular stacking. In Fig. 
[2jb) we plot the tight-binding band structure of a bilayer 
nanorribon with AA stacking. The bands also present a 
linear dispersion and they can be understood as bond- 
ing/antibonding combinations of the constituent mono- 
layer aGNR bands. 

The AB stacking can be achieved from the AA bi- 
layer geometry by displacing one graphene monolayer 
with respect to the other, in such a way that the atoms of 
one sublattice (i.e., A) of the top monolayer are placed 
over the atoms of the other sublattice (B) of the bot- 
tom monolayer. In nanoribbons, two different AB stack- 
ings are possible: 7 the AB-a stacking, shown in Fig. 
[2jc) , which yields a more symmetric geometry for infinite 
armchair nanoribbons, and the AB-/3 stacking, shown in 
[2jd). Notice that, for armchair nanoribbons, the AB- 
a configuration can be reached by displacing the top 
monolayer in the direct stacking a distance equal to the 
carbon-carbon bond acc along the ribbon length, as can 
be seen by comparing Figs. [2jb) and (c). For the AB- 
(3 stacking, the displacement is of the same magnitude 
but at 60° with the ribbon longitudinal direction, yield- 
ing a less symmetric configuration for armchair nanorib- 
bons (Fig. |2fd)). In both cases the ^4£?-stacked bilayer 
graphene nanoribbons have metallic character, and the 
conduction and valence bands have a parabolic disper- 
sion at the Dirac point. 




FIG. 2. (color online) Atomic structure geometries and 
band dispersion relations around the Dirac point for several 
armchair-terminated nanoribbons. The ribbon longitudinal 
axes are in the horizontal direction, (a) Monolayer armchair 
nanoribbon; (b) bilayer nanoribbon with AA stacking; (c) bi- 
layer ribbon with AB-a stacking; (d) bilayer nanoribbon with 
AB-/3 stacking. For this energy range, the dispersion relations 
(a)-(c) are independent of the ribbon width; case (d) corre- 
sponds to N = 17. 



B. Dirac-like Hamiltonians 

Most of the low energy properties of monolayer 
and bilayer graphene nanoribbons can be understood 
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using a k • p approxim ation, which yields a Dirac- 
like Hamiltonian! 10 * 11 * 25 ! The low-energy effective bilayer 
Hamiltonian describing the properties of a infinite AA- 
stacked bilayer has the form 



Haa — 



7i 




/ 

vpir 

7l VpTT^ 

V 71 vpir 





7i 



/ 



(3) 



where n = k, r 



iky /lC 



?k = tan -1 (k x /k y ), and 
(k x ,k y ) is the momentum relative to the Dirac 
point. The Hamiltonian acts on a four-component spinor 

(4>a\ <^b , 4$ i ^b)- The eigenfunctions of this Hamil- 
tonian are bonding and antibonding combinations of the 
isolated graphene sheet solutions, 



sfi = sv F k±j u V S A ± 



AA 



se l8k 
±1 
y ±se lB * J 



(4) 



with s = ±1. 

The low-energy Hamiltonian of the AB stacking 
reads' 



Hab = 



with eigenvalues 



^ = ^(7i± 



/ 










7i \ 




VpTT 
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7i 







/ 



(5) 



Avpk 2 + a 



s = ±1. 



(6) 



For a given eigenvalue _E,the wavefunction takes the form 



l 



\ 



E 
v F ke l9 
■^^(v F k 2 

%k 2 -E 2 
71 



71 E 



E 2 



J 



(7) 



In accordance with the geometry shown in Fig. [Tj we 
assume for nanoribbons that the system is invariant in 
the x direction, and therefore k x is a good quantum 
number. In the case of metallic aGNR, the boundary 
conditions are satisfied^ for k y —0 independently of the 
nanoribbon width; this k y —0 state is the lowest energy 
band confined in the aGNR. We have checked that the 
dispersion of the lowest energy band obtained by solving 
the Dirac model coincides with that obtained by diag- 
onalizing the tight-binding Hamiltonian for the mono- 
layer, bilayer AA and AB-a nanoribbons. Therefore, the 
Dirac approximation is a good description for the low 
energy properties of these nanoribbons, Fig. [2ja)-(c). 
This is not the situation for bilayer graphene nanorib- 
bons with AB-j3 stacking. In this case, the atomic asym- 
metry at the edges of the ribbon is not captured by the 



Dirac model, which is a long-wavelength approximation. 
Therefore, we should describe the electronic properties of 
nanoribbons with AB-/3 stacking using the tight-binding 
Hamiltonian. 



III. CONDUCTANCE 



A. Tight-Binding approach: Landauer-Biittiker 
formalism 



Due to the lack of translational invariance of the sys- 
tem, in the tight binding model we calculate the elec- 
tronic and transport prop erties using the surface Green 
function matching met hod GHHH To this end, the system is 
partitioned in three blocks: two leads, which we assume 
to be semi-infinite aGNR, and the conductor, consisting 
of the bilayer flake. The Hamiltonian is 



H = H C +H R + H L + V LC + V RC , 



(8) 



where He, Hl, and Hr are the Hamiltonians of the cen- 
tral portion, left and right leads respectively, and Vlc-, 
Vrc are the coupling matrix elements from the left L and 
right R lead to the central region C . The Green function 
of the conductor is 



Gc(E) = (E-H l 



c 



(9) 



where = Vic9tV} c is the selfenergy due to lead £ — 
L, R, and ge = (E — H/)^ 1 is the Green function of the 
semi-infinite lead fP^ 

In the linear response regime, the conductance can be 
calculated within the Landauer formalism as a function 
of the energy E. In terms of the Green function of the 



g;y 

systempHMMl it rea ds 



2e 2 

G=—T(E) 



2e 2 



(10) 



where T (E) , is the transmission function across the con- 
ductor, and = — sj] is the coupling between the 
conductor and the I — L,R lead. 



B. Continuous approximation: wavefunction 
matching 

In the low-energy limit, we can obtain the conduc- 
tance of the system by matching the eigenfunctions of the 
Dirac-like Hamiltonians. As commented above, we con- 
sider incident electrons from the lowest energy subband, 
which correspond to a transversal momentum k y =0 in 
aGNRs. Assuming an electron with energy E coming 
from the left monolayer ribbon, we compute the transmis- 
sion coefficient t to the right monolayer lead. In the cen- 
tral part the wavefunctions are linear combinations of the 
solutions of the bilayer nanoribbon Hamiltonians given in 



Sec. II B at the incoming energy E. The transmission, 
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reflection and the coefficients of the wavefunctions in the 
bilayer part are determined by imposing the appropriate 
boundary conditions at the beginning (x — 0) and at 
the end (x = L) of the bilayer region. Matching of the 
wavefunctions amounts to require their continuity. As 
the Hamiltonian is a first-order differential equation, cur- 
rent conservation is ensured automatically. The precise 
boundary condition depends both on the lead configura- 
tion (1 — >• 1 or 1 — > 2) and on the stacking. 



1. A A stacking 

In this stacking, each atom A1(B1) has an atom 
A2(B2) on top of it. The dispersion in the central part 
is given by Eq. Q, and for each incident carrier with 
momentum k x , there are always two reflected and two 
transmitted eigenfunctions with momenta ±(k x ±ji/vp)', 
see Fig. gb). 

In the 1 — > 1 (bottom-bottom) configuration the wave- 
function should be continuous in the bottom layer, i.e. 
4>^\x) continuous at x = and <j)g\x) continuous at 
x = L; for the top layer 



4 2) (x = 0) = «4 2) (x = L) = 0. 



(11) 



From these boundary conditions we obtain the transmis- 
sion 



2 AA 



sin 



4 7i L 



2 cos — COS 

Vp 



2 JiL 



a 4 2lL 

5 VF 



(12) 



In the 1 — > 2 configuration the bottom wavefunction 

(f>A (x) and the top wavefunction <j)g{x) should be con- 
tinuous at x = and x = L respectively. In addition, the 
hard-wall condition should be satisfied: 

4>f{x = 0)=4%\x = L)=Q. (13) 

The above boundary conditions yield the transmission 

,4 7i£ 



T\T = 1 - 



cos' 



Vp 



2 1 



2EL 



x 2 TlL 

Vp 



=4 ixk 

Vp 



(14) 



We see from these equations that the conductance 
changes periodically as function of the incident energy 
and length of the bilayer flake. For fixed L, the trans- 
mission is a periodic function of the incident energy. In 
the bottom-bottom geometry there are antiresonances, 
T^ 1 =0, at energies given by ^^{n + h), with n = 
0,1,2.... These energies corresponds to quasilocalized 
states in the top part of the bilayer flakeP^ The paths 
through the bottom graphene ribbon and through the 
quasilocalized state of the top flake interfere destruc- 
tively, producing the antiresonanceP*HH j n ^ ne bottom- 
top configuration, the momenta of the quasilocalized 
states of the bilayer flake are shifted in — ~, so the an- 
tiresonances occur at energies ^jf-n, with n = 0, 1, 2... . 



For fixed energy, the conductance varies periodically 
with the length of the bilayer flake. There is a period, 
nvp/E, related to the energy of the incident carrier; other 
periods are harmonics of that imposed by the interlayer 
hopping, ttvf/ji- The dependence of the conductivity 
on 71 can be understood by resorting to a simple non- 
chiral model with linear dispersion. Consider an inci- 
dent carrier from the left with momentum k x and en- 
ergy E = vpk x in the bottom sheet. When arriving at 
the bilayer central region, the incident wavefunction de- 
composes into a combination of bonding (b) and anti- 
bonding (a) states of the bilayer with momentum fc b ( a ) 
= k x ± 'ji/vp . The conductance through the bilayer 
region is proportional to the probability of finding an 
electron at the top (bottom) end of the central region, 
l±cos(fc & — k a )L = l±cos7iL/t>i?, depending of whether 
the system is in the 1 — >• 2 or in the 1 — >• 1 configuration. 
This simple model explains the dependence of the con- 
ductivity on harmonics of cos "fiL/vp and also why the 
1 — > 2 and the 1 — > 1 transmissions are in counterphase. 
The phase opposition is more evident in the E — > limit 
of Eqs. (12 1 and (14 1, which give an E = conductance 



in the bottom-top configuration equal to zero, whereas 
in the bottom-bottom configuration it has a maximum 
finite value that depends on the flake size: 



T 1 A X 1 {E = Q) = 1 



4 sin 



cos 



4 7i 
2 2 7 iL 

Vp 



T 1 A 2 2 (E = 0)=0. 



(15) 



(16) 



2. AB stacking 

In this stacking only the atoms A of layer 1 and the 
atoms B of layer 2 are directly connected by tunneling. 
The dispersion in the central part is given by Eq. ©. 
For an incident carrier with \E\ > 71 and momentum k x 
there are always two reflected and two transmitted eigen- 
functions with momentum ±fcx(2^ = ^\/k x (k x ±71/^) 
in the bilayer region, see Fig. \2\c). However, for inci- 
dent wavefunctions with \E\ < 71, there are only one 
reflected and one transmitted central wavefunctions with 
momenta ±fci = ±y // k x (k x + ji/vp)- In addition, there 
are an evanescent and a growing state with decay con- 
stants K = ±^k x (^fi/vF — k x ). Therefore, the conduc- 
tance of the system depends on whether the energy of 
the carrier is larger or smaller than the interlayer hop- 
ping. For \E\ > 71, there are two channels in the cen- 
tral region and the interference between these channels 
produces antiresonances, whereas for \E\ < 71 only an 
electronic channel is present in the central region, and 
Fabry-Perot interference can occur. Analytical, but very 
large and impractical expressions can be obtained for the 
conductance in the AB stacking. Therefore, we choose 
to present the expressions for the transmission in the low 
and high energy limit. In the next section, when com- 
paring with the tight-binding results, we plot the exact 
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results obtained from wavefunction matching in the con- 
tinuum approximation. 

The boundary conditions for AB stacking in the 
bottom-bottom configuration are similar to those of the 
AA case: (fy^'(x) and <pg\x) should be continuous at 
x = and x = L respectively, and 



<j>% ) (x = 0) = ct>% ) (x = L)=0 



(17) 



In the low energy limit, E -C 71, the AB stacking 
conductance in the bottom-bottom configuration takes 



the form 



r^ 1 (^«7i) = i-- 



4E (cos fciL+cosh kL) 2 

71 (cosh kL sin k\L— cos k\L sinh kL) 2 



which presents resonances when tanfciL = tanhreL; for 
large L this occurs when L = (n + being n an 

integer. For E — > the system has transmission unity. 

In the limit of large energy, E 3> 71 and in the bottom- 
bottom configuration the transmission is 



rAB 1 (^»7i) = l 



11 + 4cos2fciL + 4 cos (h - k 2 )L + cos2(fci - k 2 )L + 4cos2fc 2 L + 8cos(fci + k 2 )L 



(19) 



This transmission presents antiresonances associated 
with destructive interferences of the two electronic paths 
in the bilayer region. The behavior of the conductance is 
similar to that of the AA stacking, Eq. (12). There are 



periodicities associated with the energy of the incident 
electron: for E 3> 71 , 2k\L <~ 2k 2 L ~ (ki + k 2 )L ~ 
2EL/vp; and there are also periodicities associated with 
the interlayer hopping. The lower harmonic in the AB 
stacking, kl ~ k2 L ~ ^— , is half the basic harmonic in 
the AA stacking, and this reflects the fact that in the 
AB stacking only half of the atoms have direct interlayer 
tunneling. 



I 

In the bottom-top geometry <j)^2\x) and (^bHz) should 
be continuous at x — and x — L respectively, and 



4 2) 



(x = 0) = <${x = L) =0 



(20) 



In the AB stacking, interlayer tunneling connects Al 
atoms with B2 atoms; this arrangement determines the 
form of Eq. (|20|). 



For E << 71 the bottom-top transmission can be ap- 
proximated as 



7^ 2 (£«7i) = l- 



( 1 + cos fci L cosh kL + — sin k\ L sinh kL 
V 71 



4^- (cosh kL sin k^L + cos k\L sinh kL) + ( 1 + cos k\L cosh kL — 3~ sin A^L sinh 



(21) 



It can be seen that the 1 — > 2 conductance goes to zero 
when For E -> 0, T A ^ 2 tends to zero, and it is comple- 
mentary to T^g 1 . For large L, T\~g 2 presents resonances 



at L = (n + \)^~, with n integer. For energies larger 
than the interlayer hopping the conductance can be ap- 
proximated as 



7^ 2 (£»7i) = l 



> cos 



11 - 4cos2fciL - 4 cos (ki - k 2 )L + cos2(/ci - k 2 )L - 4cos2fc 2 i + 8 cos [k\ + k 2 )L 



(22) 



In this energy limit, the interference between different 
electron paths through the systems produces antireso- 
nances. Similarly to the bottom-bottom configuration, 
for a fixed energy E S> 71 the transmission varies peri- 
odically with L, with one period given by the incident 
energy, ttvf/E, and others related to the interlayer hop- 
ping, oc irv F /ji. 



IV. RESULTS 

As the systems possess electron-hole symmetry, we 
concentrate on energies E > 0. Let us recall here that 
the length of a unit cell (u.c.) for an armchair ribbon 
is 3acc = VSa-o- I n the following Figures, we choose to 
give the system length L in terms of the armchair rib- 



bon u.c. length, which is unambiguous for the discrete 
tight-binding model. Note that, in the continuum ap- 
proximation, the hard wall conditions at the edges of the 
system (x = and x = L) are set at two extra rows of 
atoms where the wavefunctions are imposed to be zero. 
This amounts to add to the system length the quantity 
a CCi which we take into account when comparing the 
continuum and the tight-binding results. 



A. AA and AB-a stackings 

As discussed in the previous Section, the expressions 
for the transmission (Eqs. @, @, @, @, (pi), and 



N = 17, E = 0.5 y 



(22|) demonstrate that the dependence with the system 
length has periodicities related to the interlayer coupling 
71 . This is evident in Fig. |3j which shows the length de- 
pendence of the conductance at a fixed energy E = 71/2, 
for the stackings AA and AB-a and the two lead config- 
urations, 1 — > 2 and 1 — > 1. Here we depict the tight- 
binding results with circles and the continuum ones with 
full lines. The tight-binding calculations are performed 
for a ribbon of width N = 17, but for this energy range 
only one channel contributes to the conductance in the 
monolayer and at most two channels in the bilayer flake, 
so the conductance is independent of N. The agreement 
between the two models is excellent for these stackings 
and energy range. As expected, the AA stacking shows 
clear antiresonances as a function of length, and the re- 
sults for the 1 — > 2 and 1 — > 1 configurations are exactly 
in counterphase. 

The results for the two configurations (1 — >■ 1 and 
1 — > 2) with AB-a stacking have an approximate com- 
plementarity; only at L — > there is zero transmission 
for the 1 — ¥ 2 case corresponding to a transmission max- 
imum for the 1 — > 1 system. The subsequent maxima 
and minima are slightly shifted, and more importantly, 
there are no zero antiresonances for finite length. As 
mentioned before, there is only one transmission channel 
in the bilayer, so although the conductance oscillates due 
to finite-size effects, there are not antiresonances for the 
AB-a at this energy. 

Fig. [4] shows the length dependence of the conduc- 
tance for another energy E — 2ji, where there are two 
conducting channels for both stackings. It is apparent the 
change for the AB-a case, which now presents antireso- 
nances with zero conductance. As to the AA stacking, 
the conductance for the 1 — > 2 configuration shows only 
one clear period of 16 u.c, whereas the 1 — ¥ 1 case shows 
also a 8 u.c. periodicity, stemming from the cos2EL/vf 
term in the conductance. The analytical expressions Eqs. 
( 12 ) and ( 14 1 allows us to verify that, for the 1 — > 2 case, 



this energy-dependent term cos47iL/vp combines with 
the other 71-dependent terms to yield a single period, 
whereas for the 1 — » 1 case, the cosAjiL/vf survives. 

Figure [5] shows the conductance G(E) as a function of 
energy for the two geometries considered and the most 
symmetric stackings, namely the AA and the AB-a, for 
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FIG. 3. (color online) Conductance as a function of the bilayer 
region length L for a ribbon of width N = 17 with AA and 
AB-a stackings, at a Fermi energy E = 7i/2. The top panel 
shows the 1 — > 2 configuration and the lower panels are for the 
1 — > 1 configuration, as schematically indicated in the upper 
left corners of the panels. The plots are labeled with the 
stackings {AA and AB-a). Red circles: tight-binding results. 
Black solid lines: continuum model calculations. 



a flake length of L = 10 u.c. (top panel) and L — 20 u.c. 
(bottom panel). The tight-binding results are depicted 
with circles and the continuum ones with lines. 

As discussed before, the most characteristic feature 
of the transmission is the appearance of Fano antires- 
onances with zero conductance. This can happen for any 
energy in the case of AA stacking because there are al- 
ways two conducting channels in the AA bilayer. On 
the contrary, for AB-a stacking, with only one channel 
for E < 71, the oscillations in the conductance are due 
to a Fabry-Perot-like effect, i.e., the interference of one 
scattering channel with itself due to the finite size of the 
structure. This is most clearly seen for the L — 20 case, 
where the AB-a stacking presents a non-zero minimum 
in the conductance in the (0,71) energy range, whereas 
the antiresonances above 71 clearly reach zero values. 
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FIG. 4. (color online) Conductance as a function of the bi- 
layer region length L for a ribbon of width N = 17 with direct 
and Bernal-a stackings, at a Fermi energy E = 271. The top 
panel shows the 1 — > 2 configuration and the lower panels are 
for the 1 — > 1 configuration, as schematically indicated in the 
upper left corners of the panels. The plots are labeled with the 
stackings (AA and AB-a). Red circles: tight-binding results. 
Black solid lines: continuum model calculations. 



FIG. 5. (color online) Conductance as a function of Fermi 
energy for a bilayer region of length L = 10 u.c. (top panel) 
and L = 20 u.c. (bottom panel). Dirac-like results: solid black 
lines correspond to the 1 — > 2 configuration and the dashed 
red lines for the 1—^1 configuration. Inside each panel, the 
top graph depicts the AA stacking, and the bottom graph 
shows the AB-a stacking data, as labeled therein. The tight- 
binding calculations are shown in colored circles. 



Notice as well the agreement with the continuum cal- 
culations when E — > 0. All the 1 — > 2 configurations 
have zero conductance in this limit. As to the 1 — > 1 
configuration, G(E = 0) has a in general nonzero value, 
which oscillates as a function of the system length, as 
described by Eq. ( 15 ) and it is seen in Fig. [5j The 
1 — > 1 and 1 — > 2 results for the AA stacking do show 
a certain complementarity: the conductance minima in 
one configuration coincide with the maxima of the other 
one. Furthermore, the periodicity of the conductance as 
a function of energy for the AA stacking due to the term 
cos 2EL/vf is evident in Fig. [5j and agrees perfectly with 
the value given by the Dirac continuum approximation, 
namely 0.16 70 for L = 10 and 0.079 70 for L = 20. 

The continuum model allows us to make a more com- 
plete characterization of the behavior of these systems. 



Fig. [6] shows contour plots of the transmission versus 
energy and bilayer flake length for the direct stacking 
and the two configurations, 1 — > 1 and 1 — > 2, given by 
Eqs. (12) and (14). We clearly see the main transmis- 



sion antiresonances with a 16 u.c. period, stemming from 
the interlayer hopping term nvp/'yi, as discussed in Sec. 
Ill B 1[ In fact, it turns out that for certain flake sizes 



L, the conductance is zero, independently of the energy. 
As this spatial period depends directly on the interlayer 
coupling strength 71 , we propose that this feature can be 
used to measure the interlayer hopping parameter, the 
value of which is still under debated by overlapping two 
nanoribbons and displacing one of them with respect to 
the other, the spatial period could be measured and thus 
71 would be obtained. 

As the variation of the conductance as a function of 
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length is so dramatic, from one quantum of conductance 
to zero, this system can function as an electromechanical 
switch, turning from the maximum conductance to zero 
by a displacement of a few A. The contour plots for the 
AA case also reveal the counterphase behavior of the 1 — > 
1 and 1 — > 2 configuration discussed previously. It is 
patent how the maxima of the conductance vs. L for 
the 1 — ¥ 2 system coincide with the 1 — >• 1 minima and 
viceversa. 
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FIG. 6. (color online) Transmission as a function of the 
energy and flake length for the AA stacking, as obtained from 
the continuum Dirac model. Top panel: 1 — > 2 configuration. 
Bottom panel: 1 — > 1 configuration. 
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FIG. 7. (color online) Transmission as a function of the en- 
ergy and flake length for the AB-a stacking, as obtained from 
the continuum Dirac model. Top panel: 1 — > 2 configuration. 
Bottom panel: 1—^1 configuration. 



agating eigenchannels in the bilayer flake. For this en- 
ergy range, the behavior is more similar to that found 
for the AA stacking, with an obvious difference on the 
spatial periods. As already mentioned in Sec. |IIIB 2[ the 
lower harmonic in the AB stacking is Sp-, thus yielding 
a longer spatial period (32 u.c.) that we attribute to the 
smaller coupling between layers for this case. 



Fig. [7] displays the contour plots for the AB-a case 
obtained within the Dirac model. There are two impor- 
tant differences with respect to the AA stacking. First, 
now there are two distinct energy regions, set by the in- 
terlayer coupling 71 . Below E = 71 , there are no antircs- 
onances because there is only one propagating channel 
at the bilayer. There are conductance oscillations, but 
not so marked as for E > ji, where the zero antireso- 
nances appear because of the coexistence of two prop- 



B. AB-P stacking 

Until this point, we have focused in the more sym- 
metric stackings, for which the continuum Dirac model 
and the tight-binding have an excellent agreement, as 
demonstrated. Now we turn our attention to the AB-(3 
stacking, which we can only model adequately with the 
tight-binding approach. This is because of the lack of 
symmetry of the ribbon edges, as it can be seen in Fig. 
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[2] (d). The atoms at the upper egde of the top layer are 
not connected to the atoms of the bottom layer, inde- 
pendently of the sublattice they belong, and viceversa. 
Such a feature cannot be well described by the contin- 
uum Dirac Hamiltonian given by Eq. |5]) , which assumes 
that all carbon atoms in the A sublattice of the bottom 
layer are connected to the B atoms on the top layer. This 
difference is not very important for wide ribbons, but it is 
noticeable for the narrow cases, for which the proportion 
of atoms at the ribbon edges is non-negligible. 

One way to assess the importance of the edge effect is 
to check the energy difference between the first and the 
second subband for E > 0. For a AB-a nanoribbon is 
always 71, whereas for AB-/3 nanoribbons it depends on 
the ribbon width, as it can be seen in Figs. [| (c) and 
(d). Size effects are related to the ratio of atoms which 
are not well described by the continuum AB Hamiltonian 
of Eq. This bring s in a dependence on the ribbon 

width, as shown in Fig. [8j depicting the conductance for 
three energies and ribbon widths TV for the two configu- 
rations, 1 — > 1 and 1 — > 2. Notice the dependence on the 
ribbon width; the conductance results demonstrate that 
size effects are still important for N s» 30. For the lowest 
energy depicted, for which there is only one propagating 
channel in the bilayer flake, the three widths show a sim- 
ilar behavior for sufficiently long flakes (L > 10). How- 
ever, for the highest energies the disagreement is patent, 
due to the dependence of the longest spatial period on 
the system width. The different periods are more clear 
for the energy E = I.571, for which at least half a wave- 
length of the oscillation can be appreciated for the three 
ribbon widths. Notice that the case E = O.571, shown in 
the central part of both panels in Fig. [HJ is also depicted 
for the AB-a stacking in Fig. [3] This striking difference 
in the conductance for the two AB stackings is due to 
the fact that in the AB-a case theres is only one channel 
for this energy, whereas in the AB-j3 there are already 
two. 
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FIG. 8. (color online) Conductance as a function of the length 
of bilayer region in AB-/3 stacking for three Fermi energies 
(E = O.271, E = O.571 and E = I.571) for three ribbon 
widths: N = 5 (dotted blue line), N = 17 (black solid line), 
and N = 29 (red dashed line). The E = O.571 and E = 
I.571 curves have been shifted up in one and two conductance 
units respectively for the sake of clarity. Top panel: 1 — > 2 
configuration. Bottom panel: 1 — > 1 configuration. 



In this work, we have studied the conductance of 
a graphene bilayer flake contacted by two monolayer 
nanoribbons. Two contact geometries have been con- 
sidered: either the left and right lead are contacted to 
the same layer of the flake or to opposite layers. Fur- 
thermore, three different stackings for the graphene flake 
have been taken into account, namely, AA, AB-a and 
AB-p. 

We have calculated the conductance with a tight- 
binding approach and also by performing a mode- 
matching calculation within the continuum Dirac model, 
by choosing the appropriate boundary conditions. We 
have explained the features in the transmission and ob- 
tained analytical expressions that allow us to elucidate 
the transport characteristics of these structures. We have 
found several periodicities on the conductance, related to 



the energy and the interlayer coupling of the system. 

In particular, for the AA configuration, we have found 
that the conductance through the flake shows Fano an- 
tiresonances, that we have related to the interference of 
two different propagating channels in the structure. For 
a flake of length L, the main transmission period is given 
by ttvf/L. For a fixed incident energy, the conductance 
as a function of the system length L oscillates with two 
main periods related to the energy E and the interlayer 
coupling 71. 

For the AB stacking, we have found two distinct be- 
haviors as a function of the incident energy E: for ener- 
gies larger than the interlayer hopping 71, the transmis- 
sions resemble those found for the AA stacking. This is 
due to the existence of two propagating channels at this 
energy range. There is, however, a difference on the main 
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period related to the interlayer hopping 71 , which is twice 
the period found for the AA stacking. This can be un- 
derstood by noticing that in the AB stacking only half of 
the atoms are connected between the two graphene lay- 
ers. For energies smaller than 71, the AB-stacked flake 
only has one eigenchannel, and the conductance shows 
resonances related to the existence of Fabry-Parot-like 
states in the system. 

The conductance of these bilayer flakes can oscillate 
between zero and the maximum conductance as a func- 
tion of length; thus, a system composed by two over- 
lapping nanoribbons can behave as an electromechani- 
cal switch. We propose that these characteristics can be 
employed to measure the interlayer hopping in bilayer 
graphene. Our results constitute a comprehensive view 
of transport through bilayer graphene flakes, clarifying 



the role of stacking, contact geometries, flake width and 
length in the conductance of these structures. 
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